Emergence of patterns in driven and in autonomous spatiotemporal systems 
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The relationship between a driven extended system and an autonomous spatiotemporal system is 
investigated in the context of coupled map lattice models. Specifically, a locally coupled map lattice 
subjected to an external drive is compared to a coupled map system with similar local couplings 
plus a global interaction. It is shown that, under some conditions, the emergent patterns in both 
systems are analogous. Based on the knowledge of the dynamical responses of the driven lattice, 
we present a method that allows the prediction of parameter values for the emergence of ordered 
spatiotemporal patterns in a class of coupled map systems having local coupling and general forms 
, of global interactions. 
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■ The phenomenon of pattern formation induced by external forcing on spatiotemporal systems, such as chemical 
£Nl ' reactions jl], |^, [j| or granular media ^, |(|, have received much attention. Similarly, there has been recent interest 
in experimental investigations of spontaneous pattern formation and emergence of collective behaviors in spatially 
extended systems of interacting dynamical elements, such as one-dimensional arrays of electrochemical oscillators 
0, chemical and hydrodynamical systems with global coupling ^ 0], and populations of chaotic electrochemical 
(— 1 ■ cells having both local and global interactions []10|. In this context, coupled map lattices have provided fruitful 
U: theoretical models for studying and predicting a variety of dynamical processes, including pattern formation, in 
spatiotemporal chaotic systems possessing different kinds of interaction topologies, such as local couplings, regular 
geometries, inhomogeneous and disordered networks, and global coupling |1 if . 

In this article we investigate the relationship between forced spatiotemporal systems and autonomous dynamical 
l—— '■ systems possessing both local and global interactions in the framework of coupled map lattices. We analyze the 
emergence of ordered patterns in forced spatiotemporal systems by using a model of a coupled map lattice subjected 
to an external drive. We show that, under some circumstances, this system is analogous to an autonomous coupled 
map system with both local and global interactions. Our approach is motivated by the observation that a globally 
• coupled map system behaves similarly to a single map subjected to an external drive and that this analogy may be 
used to describe the formation of dynamical clusters in globally coupled maps |Q . 

As a model of an autonomous spatiotemporal system, we consider the following coupled map network possessing 
both, local and global interactions 

x\ +1 = (1 - e 2 )f{x\) + % + f{x\- x ) ~ 2f{x\)} + e 2 H (x t \ xf, . . . , x?) , (1) 
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which can be compared to a one-dimensional coupled map lattice subjected to a uniform external drive, 



s\ +l = (1 - e 2 )/( S *) + -| [{f( S \ +1 ) + /(.r 1 ) - 2f{s\)] + e 2 L t . (2) 



In Eq. (|l|), x t (i) gives the state of element i (i = 1,2,..., N) at discrete time t; N is the size of the system; f(x) 
describes the (nonlinear) local dynamics; e% and e 2 are the local and global coupling parameters, respectively; and 
H(x\, x 2 , x\, . . . , x^) is a global coupling function of N variables assumed to be invariant to argument permutations; 
that is H(. ..,x\,... ,x\, ■ ■ ■) — H(. . . ,x° t , . . . , x\, . . .); Vi,j. This property of the global coupling function ensures 
that, at any time, each element of the lattice is subjected to the same influence of the global coupling term. Thus the 
interaction H provides a global feedback. As a specific example, we shall take the usual mean field global coupling 
H = jt Y]^-i f( x t)> although the analysis presented here is applicable whenever a permutable coupling function of N 
variables appears in the autonomous spatiotemporal dynamical system Eq. ([!]). In Eq. (|^), s\ is the state of element 
i in the driven lattice, /(s|) is the same local dynamics as in Eq. (|l|), €\ measures the local coupling, e 2 represents 
the coupling strength to the external forcing, and Lt is the uniform driving term. In general, Lt may be any function 
of time. We assume periodic boundary conditions and a quadratic map f(x) = 1 — rx 2 as local dynamics in both Eq. 
(|) and Eq. (§). 

The analogy between an autonomous coupled map system having both local and global interactions and a uniformly 
driven lattice arises because in the former system (Eqjl]) all the elements in the network are affected by the global 
coupling in exactly the same way at any time, and therefore the behavior of any map x\ in Eq.(Q) is equivalent to 
the behavior of any element in the driven lattice (Eq. S) with L t = H and initial conditions s z = x l . Additionally, if 



2 



the autonomous coupled map system in Eq.([y) reaches an ordered spatiotemporal pattern, its corresponding global 
coupling function H follows an ordered motion. Thus the associated driven lattice, Eq. (||), with a periodic forcing 
L t should exhibit an ordered spatial pattern similar to that of the coupled map system having both local and global 
interactions. In particular, periodic drives resulting in periodic patterns in the driven lattice Eq.(§) may be employed 
to predict the emergence of periodic patterns in autonomous systems described by Eq. (^) , regardless of the specific 
functional form of the permutable global coupling H and without doing direct simulations on the autonomous system. 

Let us analyze the dynamical response of the driven coupled map lattice subjected to periodic forcing. For a drive 
having period T, we denote the sequence of T values adopted by L t as {L\, L 2 , . . . , Lt}- For some parameter values 
and initial conditions, the resulting dynamics of the driven system (Eq.^J) may reach an ordered spatiotemporal pattern 
corresponding to lattice configurations that are periodic in both, space and time. The spatial wavelength fc and the 
temporal period p of the dynamics of the elements are defined by the relations s\ = s\ +k and s\ = s l t+p , respectively. 
In these cases it suffices to consider only the dynamics of the fc coupled elements forming a wavelength, and therefore 
the behavior of driven lattice (Eq.|J) can be analyzed by considering a reduced driven system of k coupled maps. 

An ordered spatiotemporal pattern having wavelength k and period p emerging in a periodically driven lattice can 
be characterized by the k x p matrix 



where the mth column contains the consecutive p values of the periodic response of an element m belonging to a 
spatial wavelength, (m = 1, 2, . . . , fc); and the nth row displays the values of all the fc elements in a wavelength (a 
snapshot) at the cyclic time step n, (n = 1,2, ... ,p). The driven lattice may reach different asymptotic ordered 
spatiotemporal patterns S(k,p) depending on the initial conditions; i.e. multistability is possible. 

Once a pattern S(k,p) appears in the driven lattice for some values of the parameters, it must satisfy the following 
set of fc x p nonlinear algebraic equations, 



In general, for emergent periodic patterns the ratio p/T = v is a natural number that characterizes the reso- 
nance between the period of the driving term and the resulting period of the lattice. Thus, the sequence of values 
{L\, L2, • ■ ■ , Lt} repeat themselves v times in the right hand side of Eqs.(Q). In addition, the presence of symmetries 
in a spatiotemporal pattern may reduce the number of independent variables in the above set of equations. Eqs. (Q) 
may yield several sets of solutions for the orbits <r™, however only the stable ones will be observed as asymptotic 
patterns in the driven lattice. As examples, consider the following patterns in the driven lattice, characterized by the 
given matrix and satisfying the indicated relations, corresponding to Eqs.(|]): 
(1) Frozen (p = 1) wavelength fc = 2, with constant drive L t — {Li}, 




(3) 



n+l ~ 



(1 - ei - e 2 )/(«C) + + fW 1 - 1 )} + *»L, 



(4) 



S(2, 1) = (a b) ■ 



(5) 



0= (1 

b= (1 



ei - ea)/(a) + ei/(6) + e 2 Li 
e 1 -e 2 )/(6) + e 1 /(a) + e 2 L 1 . 



(6) 



(2) Out of phase, wavelength k — 2 and period 



p = 2, with constant drive {£1}, 




a= (1 
b= (1 



ei-e 2 )f{b) + e 1 f(a) + e 2 L 1 
ei - £2)/ (a) + ei/(&) + £2^1 • 



(8) 



(3) Wavelength k = 2 and period p = 4, with period-two drive {Li, L 2 } 




(9) 



3 



c = (l-e 1 -e 2 )f(a)+e 1 f(b) + e 2 L 1 

d= (l-e 1 -e 2 )/(6) + e 1 /(a)+e 2 L 1 

a= {l-e 1 -e 2 )f{d)+e 1 f(c) + e 2 L 2 

b= (l-ei-e 2 )/(c)+ei/(d)+e 2 i 2 . 



(10) 



(4) Wavelength k = 3 and period p = 3, with constant drive 



5(3,3) 



( 



a b c 
b c a 
cab 



) 



(11) 



a= (1 
6= (1 
c= (1 



ei-e 2 )/(c) + f [/(6) + /(a)] + e 2 i 1 
ei -e 2 )/(a) + f [/(c) + /(&)] +e 2 Lx 
ex - e 2 )/(6) + f [/(a) + /(c)] + e 2 Lx . 



(12) 



For the above examples, the orbits a, 6, c, d can be obtained as functions of the parameters ei, e 2 , r, and the values 
Li and L 2 . 

In practice, if we are searching for the orbits <r™ in a specific stable pattern S(k,p), we can construct a bifurcation 
diagram of a periodically driven lattice of N coupled maps (N multiple of k) as a function of some parameter of the 
system, and look for windows where that pattern appears in the diagram. This may require exploring for appropriate 
initial conditions. Figures l(a)-(d) show bifurcation diagrams of driven lattices as a function of either ei or e 2 , 
where the patterns given in examples (l)-(4), additionally to other spatiotemporal states, arise. Fig. 1(a) shows the 
bifurcation diagram of a lattice subjected to a constant drive L\ — 0.39 as a function of the coupling parameter of the 
drive e 2 , with fixed local coupling ei = 0.54. As e 2 is varied, the steadily driven system exhibits several spatiotemporal 
states, such as synchronization (both chaotic and periodic), a pattern S*(2,l), and quasiperiodic k — 2 wavelength. 
The stationary orbits a and b (for either the odd or even elements of the lattice) in the pattern S(2, 1) arise in the 
bifurcation diagram of Fig. 1(a) as functions of the parameter e 2 . The values of the instantaneous mean field of the 
wavelength in each pattern S(k,p), given by (/}„ = i Yl m =i are a ^ so plotted on the corresponding regions of 

the bifurcation diagrams. In the case of Fig. 1(a), we get (/) = [/(a) + f(b)]/2, independent of time, on the region 
5(2, 1). Note that at e 2 = 0.56, the mean field (/) is equal to the value of the drive L\ = 0.39 in Fig. 1(a). Since 
(/) also corresponds to the mean field of the entire lattice, it may be inferred that for the values of the parameters 
r = 2, ei = 0.54, and e 2 = 0.56, an autonomous coupled system described by Eq. (|l|) with mean field global coupling 
may also exhibit a spatiotemporal frozen pattern with k = 2 and p = I, while sustaining a constant value of its 
global coupling function at H = 0.39. For those parameter values, the driven lattice Eq. (|J) and the autonomous 
system Eq. ([!]) are equivalent. Note that this prediction is being made without direct numerical simulation on the 
autonomous spatiotemporal system. Moreover, the intersection of (/) with the constant value of L\ in the diagram 
of the driven lattice in Fig. 1(a) readily provides a set of initial conditions for the odd and even elements x\ in Eq. (|l|) 
for observing the spatiotemporal pattern S(2, 1) in the autonomous coupled system at the parameter values e\ = 054, 
e 2 = 0.56, and r = 2; that is, x\ = a = 0.017, for i even; and x\ — b — 0.776, for i odd. Similar predictions of 
parameter values and initial conditions for the emergence of the spatiotemporal patterns S(2, 2), S(2, 4), and ,9(3, 3) 
of examples (2)-(4) in the autonomous system Eq.lQ) can be made from the intersections of the curves (/)„ with 
the values of L t in the diagrams of the periodically driven lattices of Fig. 1(b), 1(c) and 1(d), respectively. Figs. 
2(a)-2(d) display the spatiotemporal patterns S(k,p) of examples (l)-(4) emerging in the autonomous system Eq.(Q) 
while sustaining either a constant or a periodic mean field coupling function, at parameter values and orbits predicted 
from the analogy with the driven coupled map lattice. 

These results suggest that a complete equivalence between a lattice driven with period T (Eq.|^), that gives rise to 
a spatiotemporal pattern S(k,p) (p/T = v), and an autonomous system (Eq.|l|) can be established when the following 
conditions are satisfied 



where n = 1, . . . , T; and where the arguments of H are the k elements in the nth row of the matrix S(k,p), repeated 



Eqs. (y_3|) predict that the autonomous system (Eq.|lj) with local and global interactions possesses a state characterized 
by the spatiotemporal pattern S(k,p) with orbits x™ = a™{L^) 1 and by the periodic global coupling function motion 



L n — H{. . . , a x n 




(13) 



N/k times 
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H = L\,H = Ll,...,H = L* T . As an illustration, consider the pattern 5(2,1) whose orbits a(r, e\, 62, £1) and 
b(r, ei, £2j Li) can be obtained from Eqs. (0). Then Eqs. (|u]) reduce to just one equation, which in the case of a mean 
field coupling function H — J2i f{ x t)i yields the solution 

L = - arfa+k-D - (14) 

Thus for the set of parameters satisfying the above relation, the pattern 5(2, 1) can emerge in both, a steadily driven 
lattice with L t = {L|} and an autonomous system sustaining a constant mean field coupling H = L\. 

Fig. 3 displays the function H(a™) corresponding to mean field coupling as a function of the constant drive L\ for 
the patterns of examples (1), (2) and (4). The intersections of the curves with the diagonal give the solutions L\ to 
Eq. (|l3| ) for each pattern indicated. Note that Eqs. ( [l3|) can be used to predict if a given pattern S(k,p) observed 
in the driven lattice may emerge in diverse autonomous systems possessing different functional forms of the global 
coupling function H . For instance, one may ask if the pattern 5(2, 2) of example (2) can also exist in an autonomous 

system Eq. (|l|) with a constant value of a global coupling function given by the geometric mean, H = Jjf la^l 1 W . 
Fig. 3 shows H(a™) for this form of global coupling associated to the pattern 5(2, 2) (slash-dotted curve) as a function 
of Li, giving a solution L* = H = 0.3418 at the intersection with the diagonal. This prediction has been verified for 
this autonomous system. 

In summary, we have studied the emergence of induced ordered patterns in forced spatiotemporal systems by using 
a model of a one-dimensional coupled map lattice subjected to an external drive. Under some circumstances, this 
system is analogous to an autonomous coupled map system possessing a similar local coupling and an additional global 
interaction that acts as a global feedback. Thus, by exploring the dynamical responses of a driven spatiotemporal 
system, one can get insight into the conditions for the emergence of specific patterns in a class of autonomous 
spatiotemporal systems. Once an ordered pattern S(k,p) appears in a driven lattice for some values of the parameters, 
the same pattern is expected to spontaneously arise in a family of related autonomous coupled map systems that 
satisfy Eqs. jl^), which constitute the link between both systems. Conversely, a pattern observed in an autonomous 
system with both local and global interactions should also exist in an associated driven system with similar local 
couplings for some appropriate period of the external drive. This equivalence could be tested in experimental settings 
of interacting dynamical elements [0, ||, ^, (To). The method has been applied to some simple cases; however it can 
be used for more complicated patterns. Although we have considered one-dimensional diffusive local couplings, the 
analogy between a uniform external drive and a global interaction can be applied to any network of coupled maps. The 
relation between extended systems subjected to non-uniform driving fields and spatiotemporal autonomous systems 
is an interesting problem for future research. 

This work was supported by Consejo de Desarrollo Cienti'fico, Humamstico y Tecnologico of Universidad de Los 
Andes, Merida, Venezuela. 
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FIG. 1: Bifurcation diagrams of the orbits s\ of the driven lattice (Eq.^) with size N = 30 and fixed local parameter r = 2. 
The values of the periodic drive Lt are shown with dashed lines. Orbits a, b, c, d, corresponding to patterns S(k,p) of examples 
(l)-(4) are indicated. The mean (/)„ is drawn with thick lines on each region where a pattern S(k,p) appears. Regions where 
synchronization occurs are identified; regions of quasiperiodic behavior are labelled QP, and those of spatiotemporal chaos are 
labelled STC. (a) e\ — 0.54, and constant drive L\ = 0.39; bifurcation parameter is e-x. (b) t\ = 0.05, L\ = 0.36; bifurcation 
parameter 62. (c) ei = 0.54, and period-two drive Li = 0.77,1/2 = 0.22; bifurcation parameter 62. (d) 62 = 0.065; Li — 0.08; 
bifurcation parameter is e\. 
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FIG. 2: Spatiotemporal patterns on a grey scale in the autonomous system (Eqjl|) with global mean field interaction H, at 
parameter values predicted from the analogy with the driven lattice in Fig. (). Size N = 30; spatial index of x\ runs horizontally 
and time runs from bottom to top. (a) Pattern S(2, 1); r — 2, ei = 0.54, £2 = 0.56; constant H — 0.39. (b) Pattern 5(2,2); 
r = 2, ei = 0.05, £2 = 0.2; constant H — 0.36. (c) Pattern S(2, 4); r = 2, ei = 0.54, £2 = 0.51; Ji" oscillates periodically between 
the values 0.77 and 0.22. (d) Pattern S(3, 3); r = 2, ei = 0.001, e 2 = 0.065; constant # = 0.08. 
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FIG. 3: Global coupling functions H(a™) (r.h.s of Eq. (|l3|)) associated to patterns S(k,p) with constant drive, as a function of 
Li . Continuous curves correspond to mean field global coupling with the parameters of Fig. (2) in each case. Dash-dotted curve 
displays a geometric mean global coupling function associated to the pattern S(2, 2) with parameters r = 2, t\ — 0.005, ei = 0.4; 
the intersection occurs at L\ — 0.3418. 



